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Abstract 

Using  the  recently  developed  method  of  recurrence  relations 
we  have  studied  the  time  and  frequency  dependent  behavior  of  the 
two  dimensional  electronic  system  at  long  wavelengths.  Several 
families  of  autocorrelation  functions  including  the  relaxation 
and  memory  functions  are  obtained  exactly.  In  addition,  other 
linear  response  quantities  e.g.  the  density-density  response 
function  and  dynamic  structure  factor  are  obtained.  This  work, 
contains  a  detailed  application  of  the  method  of  recurrence 
relations. 


PACS:  72.10  Bg,  72.20.-i,  73.40QV 


I .  Introduction 


The  electron  gas  has  had  a  long  history  as  a  model  for 

metals  and  has  attracted  considerable  theoretical  efforts  over 

the  years. ^  Interest  in  this  model  has  been  substantially 

heightened  recently  by  the  physical  realization  of  an  essentially 

two  dimensional  electronic  system  in  MOS,  heterostructures  and 

superlattices.  '  In  addition,  the  realization  that  many-body 

effects  are  in  some  ways  more  pronounced  in  two  dimensions  than 

in  three  dimensions  has  further  given  impetus  especially  to 

theoretical  efforts.  The  static  properties  of  the  two  dimensional 

version  of  the  electron  gas  model  have  already  been  extensively 
4 

studied.  The  correlation  energy,  for  example,  has  been  calcu¬ 
lated  by  techniques  which  are  basically  similar  to  those  employed 

5-10 

for  the  three  dimensional  case.  The  dynamic  properties  have 

also  received  some  attention,  but  largely  limited  to  the  frequency 

dependent  density  response  function  calculated  by  RPA  and  other 

11-14 

perturbative  techniques.  To  our  knowledge,  there  are  as  yet 

no  measurements  of  the  dynamic  structure  factor  with  which  these 
calculations  can  be  compared. 

In  the  two  dimensional  electron  gas  model,  one  has  an 
attractive  possibility  of  finding  exact  time  and  frequency 
dependent  solutions,  normally  precluded  in  the  three  dimensional 
model.  Such  solutions  can  contribute  to  our  understanding  of 
nonlinear  behavior  of  fluids,  e.g.,  memory  effects  and  long 
time  tails  in  the  velocity  autocorrelation  function.  These 
studies  are  attracting  considerable  current  activities . ^ 

The  time  evolution  of  a  dynamical  variable  e.g.,  density- 


fluctuation  operator  may  be  studied  by  the  method  of  recurrence 
16  17 

relations.  '  Given  certain  static  properties  known  as  the 

recurrants ,  one  can  obtain  the  exact  time-dependent  behavior  of 

the  density  fluctuations,  from  which  the  associated  relaxation 

18 

and  memory  functions  as  well  as  the  random  force.  This  method 

has  already  been  successfully  applied  to  several  formal  and 

19 

physical  models.  We  present  here  a  complete  account  of  the 

time  and  frequency  dependent  behavior  of  the  two  dimensional 

electron  gas  at  long  wavelengths.  A  preliminary  account  of  this 

20 

work  was  previously  reported. 


Time  Evolution  of  Density  Fluctuations  in  the  Electron  Gas  at  T=0 


We  consider  the  two  dimensional  electron  gas  model  defined 
by  the  Hamiltonian  H 

H  ‘l  -y.Ck+ck  *  ''l  VkV*+Cg-k+CpCq  =  Ho  *  V  (1) 

k  k,p,q^ 

where  c^  &  c^  are,  respectively,  the  creation  and  annihilation 

2  2 

operators  at  wave  vector  k,  =  k  /2m,  v^  =  2-rre  /k,  where  m  and  e 
are,  respectively,  the  electron  mass  and  charge,  ft  =  1.  Our 
system  is  imposed  under  an  external  perturbing  potential  of  the 
following  form: 

Hext  *  I  ck(t)  eiUt  <2> 

k 

where  pk  is  the  Fourier  component  of  the  external  electric  field 
sufficiently  small  as  to  permit  the  use  of  linear  response  theory, 
w  is  the  frequency  of  the  field,  pk  is  the  density  fluctuation 
operator  defined  as 


P>  =  l 


Pj^  (  l  ®  Pk  6  ( 3b) 

with  p^CO)  =  Pk*  Given  an  explicit  knowledge  of  time  dependence 
in  pk(t) ,  one  can  proceed  to  calculate  e.g.  relaxation  function 
(Pk(t),  Pk)/(pk,Pk) ,  from  which  other  physical  quantities  via 
linear-response  relations.  Here  the  inner  product  means  the 
Kubo  scalar  product  (see  sec.  Ill)  . 

It  is  important  to  note  that  the  calculations  of  time 
evolution  are  carried  out  with  the  static  properties  e.g.  k  fixed 
at  some  given  value.  In  this  work,  we  shall  confine  ourselves 
to  k/k  <<  1,  where  k_  is  the  Fermi  wave  vector,  such  that  the 

F  r 

density  fluctuations  are  those  which  occur  very  near  the  Fermi 


surface.  Physically  it  means  that  excitations  in  (1)  are  limited 
to  electron-hole  pairs  only  (i.e.,  no  multipair  excitations 
allowed) .  In  this  static  regime  (k/k^  <<  1  and  T  =  0)  the  static 
susceptibility  (p^,  Pk)  and  other  static  quantities  of  interest 
are  known  exactly. 


III.  Method  of  Recurrence  Relations 


If  the  electron  gas  is  momentarily  perturbed  by  an  external 
field  coupled  to  density  fluctuations  (2)  in  accordance  with  linear 
response  theory,  the  system  will  undergo  a  relaxation  process. 

The  time  evolution  of  density  fluctuations  satisfies  the 
Heisenberg  equation  of  motion  viz.. 


Pk(t)  =  i  [H,  pk(t)  ]  . 


(4) 


One  can  put  it  in  the 

t 

dpk(t)/dt  +  /dt' 


18 

form  of  the  generalized  Langevin  equation 
Mk(t-f)  Pk(t‘)  =  Fk(t)  (5) 


where  M  and  F  are,  respectively  the  memory  function  and  generalized 
random  force.  This  form  has  the  advantage  of  introducing  the  two 
physical  quantities.  The  formal  solution  for  the  time  evolution 
of  is  given  by  the  method  of  recurrence  relations1®. 

According  to  this  method,  one  constructs  S  a  d-dimensional 
Hilbert  space  of  pk  spanned  by  a  set  of  orthogonal  basis  vectors 
fQ  f^...f(,_^.  This  physical  space  S  is  an  abstract  space  realized 
by  the  Kubc  scalar  product 

(X,Y)  =  JSdu  <euH  X  e“uH  Y+>  -<X><Y  +  >  (6) 


for  every  X,  YCS,  where  <...>  denotes  the  canonical  ensemble 
average, +  Hermitian  conjugation  and  8  is  the  inverse  temperature. 
The  method  of  recurrence  relations  shows  that  in  such  a  space 
the  basis  vectors  satisfy  the  following  recurrence  relations  (RRI) : 


fv+l  =  £v  +  fv-l'  0  i  v  d_l 


where  fy  =  i  [H,  f y ] 


(7) 

(8) 


and 


=  <Vfv)/{fv-l'  fv-l)- 


(9) 


Here  A  ,  a  relative  norm,  will  be  referred  to  as  the  vth  recurrant . 


The  boundary  conditions  require  that  =0  fi  a  =  1. 

By  choosing  f  appropriately,  one  can  readily  obtain  all 

other  fv's  by  the  RRI.  Given  these  basis  vectors  for  S,  the 

solution  of  the  Heisenberg  equation  of  motion  for  p^(t)  may 

be  formally  given  an  orthogonal  expansion  as 

p  <t)  =  j  av<t)f  UO) 

v"  0 

where  ay(t) 's  are  real  time-dependent  functions  i.e.,  auto¬ 
correlation  functions.  If  in  fact  one  chooses  fQ  =p^/  a^ft) 
corresponds  to  the  relaxation  function  of  linear  response  theory. 
The  most  important  property  for  these  functions  is  that  they 
satisfy  a  lecurrence  relation  (RRII) : 

Av+1  av+l(t)=  "  av(t)  +  av-l(t)'  0  -  v  -  d_1  (11) 

where  a  (t)  =  da  (t)/dt  and  a  ,  (t)  =  0.  With  the  choice  f  =  p. 

the  boundary  conditions  are  that  ao(t=0)  =  1,  ay(t=0)  =  O  if  v>  1. 
The  above  recurrence  relation  is  realized  by  the  recurrants  A^'s, 
which  are  static  properties.  For  the  two  dimensional  gas  at  long 
wavelengths  (k/kF <<  1)  and  T  =  0  they  are  exactly  calculable. 

Hence,  it  may  be  possible  to  solve  the  RRII  to  obtain  the  auto¬ 
correlation  functions  a^ ( t)  . 

The  random  force  for  is  a  vector  in  a  subspace  of  S . 

It  may  also  be  given  an  orthogonal  expansion 
d-1 

Fk(t)  -  lMl  fv  (12) 

where  b^(t) 's  are  real  time-dependent  functions,  of  which  b^(t) 
is  related  to  the  memory  function  (t)  =  b^(t).  These  auto¬ 

correlation  functions  also  satisfy  a  recurrence  relation  like  RRII 


(13) 


AV+1  bv+l(t)=  ■  bv(t)  +  Vl(t)'  1  -v-  d_1* 

Observe  that  does  not  enter  into  (13).  Hence,  b^ttj's  do 
not  depend  on  the  first  recurrant.  The  boundary  conditions 
require  that  b^(t=0)  =  1  &  bv(t=0)  =  0  if  v  >  2.  The  two 
families  of  autocorrelation  functions  are  related  by  a  convolution 
relation 

ay(t)  =  convol.  (bv  ( t)  *  ao(t)),  l<v<d-l.  (14) 

Hence,  if  aQ(t)  is  known,  ay(t)  is  knowable  by  the  RRII  and 
b^(t)  via  convolution.  The  method  of  recurrence  relations  is 
exact.  If  the  recurrants  are  exactly  given  in  some  static  regime, 
the  RRII  yields  exact  time  dependence  in  that  static  regime. 

For  example,  the  response  function,  one  of  the  autocorrelations 
av(t)'s,  satisfies  the  moment  sum  rules  exactly  to  all  orders 
in  the  static  regime  where  the  recurrants  are  valid. 


IV.  Basis  Vectors  and  Norms 


fo  =  pk 
fl  =  pk 


By  choosing  f  =  p,  we  can  directly  obtain  from  the  RRI  all 
other  basis  vectors  spanning  S: 

(15a) 
(15b) 
(15c) 
(15d) 

Pk  (15e) 

1  A4  + 

etc,  where  the  number  of  dots  or  simply  the  number  above  pk  refers 
to  the  number  of  times  the  commutator  of  pk  with  H  is  being  nested 


f  2  = 

Pk 

+ 

41  pk 

f  3 

o'. 

■  X 

X 

( A^  + 

A  \  A 

^2  uk 

f  4  - 

X 

k 

+ 

+ 

r—i 

<1 

A2  +  L'^  *p’k  +A1  A3 

f  5  " 

6. 

+ 

(Ai  + 

A2  +  A3  +  A4^k  +  ^A 

e.g 

•  •  • 

*  '  pk 

—  X 

[H 

'  PkJ 

=  i3 

[H 

,  [H,  (H,  Pk] 

]  } 

=  a 

k‘ 

The 

norms 

of 

these  vector 

s  are : 

<V 

V 

= 

tpk' 

Pk)  = 

xk 

<V 

V 

= 

V 

(f2. 

V 

= 

(pk' 

pk}  “ 

(A1 

<»*-  pk> 

(f3. 

f3} 

= 

{pk' 

pk} 

(A1 

+  &2)  (p'k,p 

k> 

<f4' 

f4> 

= 

(^k' 

V  - 

(A1 

+  a2  +  ^3) 

(Pk, 

pk> 

+  A1A3 

( f  5 ' 

f5> 

= 

<*k' 

A  )  - 
pk' 

<A1 

+  A2  +  A3  + 

A4> 

'  V 

+  (A1A3+A1A4 

+  i244) 

r«vpk> 

etc 

.  ,  where 

we 

have 

used 

the 

identity  (p. 

mv 

P)  = 

-  -1? 

-1  m+l. 
P  ,  P  ) 

(16a) 

(16b) 

(16c) 

(16d) 

pk>  < 

<16f> 


,n+l  m  ■* 
(  P  / 


n  or  m  >  1 . 


Observe  that  to  calculate  f  or  its  norm,  we  need  to  know 

the  norms  of  f  . ,  f  _,  ...f  .  Thus,  if  we  calculate  the  norms 
v-1  v-2  o 

successively  starting  with  the  basal  vector  fo,  we  can  obtain 
all  the  basis  vectors  and  their  norms  explicitly.  For  the 
electron  gas  at  long  wavelengths,  these  norms  can  be  obtained 


i  a 


by  working  out  the  nested  commutators  of  p^  and  evaluating  the 
inner  products  at  T  =  0.  Observe  also  that  =  i[H,pkJ  =  i  [Hq, 
p^]  =  p^°^  /  where  Hq  is  the  noninteracting  part  of  H.  As  a 
result,  the  noninteracting  or  ideal  Hamiltonian  will  have  an 
important  input  in  the  behavior  of  the  recurrants. 

Ideal  system 

Here  we  express  k  in  units  of  kp  such  that  Gp  =  l/2m 
where  £p  is  the  Fermi  energy  in  two  dimensions.  In  these  units, 
the  norms  for  the  two  dimensional  electron  gas  are  given  below: 


(fo'fo>  <0)'  Xk10'-  k-<2 

(17a) 

(f^fj.)  r°}=  k2  Gp/TT 

(17b) 

(f2,f2)  (o)  =  (k4  +  k6)  gf2A 

(17c) 

(f3,f3)  (0)  =  (k6  +  4k8  +  0  (k10)  )  Gp5/ *T 

(17d) 

(f4,f4)  (o)  =  (k8  +  10k10  +  0  (k12)  )  gf7/tt 

(17e) 

(f5,f5)  (o)  =  (k10  +20k12  +  0(k14))  ef 9/tt 

( 1 7  f ) 

etc.  For  the  recurrants  valid  in  the  static  regime  of  k«  1,  we 


shall  retain  only  the  leading  order  in  k: 
,  (o) 


2  2 
-  2k  eF 


(18a) 
(18b) 

In  the  two  dimensional  ideal  electron  gas  at  long  wave- 


(o) 


A  (0)_ 

- 


,2  2  _  . 

k  Gp  = 


lengths  the  norms  are  finite  and  non-vanishing  for  all  v’s.  Hence 
the  dimensionality  of  the  Hilbert  space  of  p^  is  infinite  (d=°°)  . 
Interacting  systems 

For  the  static  regime  of  small  k,  we  find  that  the  norms  are 


to  the  leading  order  in  k  as  follows: 


(f  ,f  )  =  >,  (o)/(l  +  v.  >:v(o)) 


(19a) 


Hence,  to  the  leading  order  in  k, 

A1  =  2A  +  r  (20a) 

Ay  =  Ay(o) =  A,  v  >  1  (20b) 

where  f  =  (u  x)  =  2ttp  e  k/m,  where  p  is  the  number  density. 
For  the  two  dimensional  interacting  electron  gas  at  long  wave¬ 
lengths  the  dimensionality  of  the  Hilbert  space  of  is  also 
infinite. 


- 

>vw\v 

■  ■  k  *  I 


.  * .  v 

-  y  '  W  - 


y"*  «>  y> 

V.n'  i 
V\'- 


Autocorrelation  Functions  for  Ideal  System 


Using  the  recurrants (18a, b)  obtained  for  the  ideal 

electron  gas,  we  can  now  obtain  the  autocorrelation  functions 

av^°^  (t)  &  b  ^  (t)  from  the  RRI  valid  at  long  wavelengths. 

2 

Defining  A  =  k  y  ,  we  can  decompose  the  RRII  as  follows: 

h  V2  (t)  =  -  aQ(o)  (t)  (21a) 

*  y2  av+1(0) (t)  =  -  a^o) (t).  +  av_1(o) (t),  v  >  1.  (21b) 

We  look  for  a  function  which  satisfies  the  above  equations.  Let 
ay(o) (t)  -  X~v  Uv  hv  (Xt)  (22) 

where  X  is  a  scale  factor,  Uv  is  a  coefficient,  and  h^  is  a 
function,  all  to  be  determined.  The  boundary  condition  a  ^ (t=0) 
imposes  that  Uq  =  1  &  hQ(0)  -  1.  Substituting  (22)  into  (21a, b) 
we  get 

h  y2  X~2  hx  =  -  hQ»  (23a) 

'  w2  hv+l  -  -  <Wl>  V  +  Vl’  v  i  <-23b> 

where  h'  =  ch(x)/dx.  Clearly  the  scale  factor  is  X  =  y,  which 
somewhat  simplifies  (23a, b)  to 


h  U.  h.  =  -  h  '  (24a) 

11  o 

‘V^Vl’  hv+l  ”  -  <Wl>  V  +  Vl'  v  i  x-  (24b) 

Now  the  above  may  be  compared  with  the  recurrence  relation  for 

21 

the  Bessel  function 


J,  =  -  J  ’ 
1  o 


J  ,  =  -  2  J  1  +  J  ,  , 
v+1  v  v-1 ' 


v  >  1  . 


(25a) 

(25b) 


If  U  =  2V,  (24a, b)  become  identically  reduced  to  (25a, b) .  Hence 


av(o)  (t)  =  2V  y_v  Jv  (yt)  ,  v  >0. 


(26) 


Observe  that  the  solution  satisfies  the  boundary  condition 


that  a. 


(o) 


(t=0)  =  0,  v  >  1.  Thus  the  RRII  realized  by  the 


recurrants  of  the  two  dimensional  ideal  electron  gas  is  isomorphic 
to  the  recurrence  relation  for  the  Bessel  functions.  The  short 
and  long  time  behavior  of  the  relaxation  functions  ay ^  (t)  is 
given  directly  by  the  asymptotic  properties  of  the  Bessel 
functions  e.g. 

„  (o) 

o  (t  -*>»)'•  t  cos(yt  -  *jtt)  .  (27) 

For  dynamical  analysis  it  will  be  useful  to  have  the 

Laplace  transform  of  the  relaxation  functions: 

av(o) (z)  =  T  (av(o) (t) ]  (28) 

where  T  is  the  Laplace  transform  operator. 

The  autocorrelation  functions  b  ^  (t)  can  be  obtained  by 

the  convolution  eq.  (14)  which  when  Laplace-transformed  takes 


a  simple  form: 

bv(z)  =  av(z)/aQ'(z)  (29) 

where  b^(z)  =  T  [b^  ( t)  ]  ,  v  >  1.  Using  (28)  we  obtain  directly, 
bv(o)(z)  =  2V  y~2  ((z2  +  y2)4  -z)V,  v  >  1.  (30) 

Hence, 

bv(o)  (t)  =  T'1  [bv{o)  (z)  3 
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=2V  y1  V  v  J  (yt/yt)  ,  v  >  1.  (31) 

v 

The  memory  function  for  the  two  dimensional  ideal  electron  gas 
becomes 

M(o)  (t)  =  A1(o)b1(o)  (t)  =  y2  J1  (yt/yt)  (32) 

~ t  ^^2  cos  ( y t  -  3t/4)  if  t  -*■  00 .  (32a  ) 


Other  linear  response  quantities  e.g.,  the  response  function 
Xjc(w),  the  dynamic  structure  factor  S^to)  can  now  be  directly 
obtained  from  the  autocorrelation  functions.  From  linear 
response  theory  (henceforth  suppressing  k  dependence) 


(34) 


Hence , 

X  ^  ( t)  =  y  J1  (yt) 

and 

X  (o,(«)5T  [X(o)(t)]2=iu 

2  2  -k 

=  l-iw(u  -  cj)  ,  0  <  w  <  y  (35a 

=  1  -  u(cu2  -y2)  **  ,  y  <  a)  <  ».  (35b 

This  is  equivalent  to  the  well  known  result  due  to  Stern1^  at 
long  wavelengths.  Also  since  it  S(u)  =  -  Ira  \(.z  =iaj)  ,  we  obtain 
directly 

S(o)  (a)  =  S(o)  (a)/x(0) 

=  co(y2  -  o2r^  0  <  u  <  v  (36a 

-  O  ,  y  <  a  <  “>f  (36b 

Observe  that  (u)  has  a  power  law  singularity  at  u  » y- , 

suggestive  of  an  apparent  collective  excitation.  The  dynamic 
structure  factor  satisfies  the  standard  moment  sum  rules,  to  all 
orders  at  long  wavelengths.  In  this  static  regime  the  time 
evolution  in  the  two  dimensional  ideal  electron  gas  is  thus 
exactly  and  completely  characterized. 


VI.  Autocorrelation  Function  for  Interacting  System 

It  was  shown  that  for  the  two  dimensional  interacting 
electron  gas  at  long  wavelengths  Ay=  v  >  2,  although 

A1  ^  (see  eq.  20a,b) .  Since  the  autocorrelation 


functions  b^,(t)  depend  not  on  the  first  recurrants  but  on 

others  only  we  must  have  b(J(t)  =  b  ^  (t)  ,  v  >  1.  We  take 

advantage  of  this  simplification  to  obtain  the  autocorrelation 

functions  ay(t)  as  follows:  By  taking  the  Laplace  transform 

of  the  RRII,  we  get  for  v  =  1, 

A  b.  (z)  =  (a  ( z) ) -1  -  z.  (37) 

v  1  o 

Using  (30)  for  b.^(z)  =  b^°^  (z) ,  we  get 


aQ(z)/sy  = 

(z2  +  u2)*5  -  (1  -  s)z 

(38) 

1  2  ^  2 
y  +  az 

where  s  =  A.^0^ 

/Ax  =  y2/  (u2  +  2D 

(39) 

a  e  l  -  (l 

-  s)2  =  (x2  +  D/(x2  +  >j)2 

(40) 

where  x  =  f/u  . 

Note  that  0  <s  <_1  and  hence  0  <a  <1. 

In  both 

cases,  the  upper  bound  is  reached  in  the  ideal  limit.  By  taking 


an  inverse  transform  we  have 


a  (t) 
o 


2iri  J 


dz  e 


zt 


2  2 

(z  +  y  )  ’  -(1  -  s)z 


+  az 


(41) 


where  the  contour  c  runs  along  the  imaginary  axis.  There  are 
two  types  of  singularities  on  the  imaginary  axis.  There  is  a  cut 
from  z  =  -  iy  to  z  =  iy.  In  addition,  there  are  a  pair  of 
isolated  poles  beyond  the  cut  at  z  =  ±  ia”*1  y  =  ±  iw^.  *n 


Appendix  A  we  have  evaluated  this  integral: 


ao(t)  =  A s  l  (-a)n  (3/3ytrn  [J^y^/yt]  +  A  cos  to t  (42) 


where 


A  =  1  -  (1  -  a)’ 

5 

Ap  ■  2  [  (1  -  a) %  -  (1  -  a)  3 /a. 


(42a) 

(42b) 


We  can  now  check  various  boundary  conditions  on  aQ(t)  . 

In  Appendix  B,  we  show  that  aQ(t=0)  =  1 .  In  Appendix  C,  we 

show  aQ(t)  =  JQ(yt)  at  the  ideal  limit  a=  1.  It  is  easy  to  show 

that  as  a-*-  0,  a  (t)  =  (1  -  >*a)  cos  u  t  +  >ja  J.  (yt)/yt  +  O  (a2). 

o  pi 

For  t  a  (t)  -  t-3//2  cos  (yt  -  3tt/4)  +  A  cos  to  t.  Also  given 

o  P  P 

a  (t)  ,  the  other  a  (t)  's  can  be  obtained  from  it  with  successive 
o  v 

applications  of  the  RRII,  e.g.,  Aj^a^Ct)  =  -  ao(t),  A2a2(t)  =  -  a^ 
+  aQ(t) ,  etc. 

In  Appendix  D,  we  have  obtained  a  closed  form  expression 

for  the  Fourier  transform  of  a  (t)  using  (42) : 

o 


a  (u)  = 
o 


As  (y2  -  J)h 

-rr(y  -  aw  ) 


'  0  <  u)  <  u 


(43a) 


=  “s  A  [6(u  -  u  )  +  6(ui-+  (!)  )  ]  ,  y  <  <0  <  “.  (43b) 

P  P  P 

In  Appendix  E,  we  have  obtained  the  response  function 
Xk(w)  =  (w) /Xj,  using  (42):  Again  suppressing  k  dependence. 


As  (1  -  a) h  u2 

Re  x  («*>)  =  1  +  - 2 - 2 — 

(y  -  Guo  ) 


>  0  <  to  <  y 


„  i.  +  As  i <  1  -  a ) 11  (i  -  v2/^2) ‘’i^2  ,  v  <  „  <  a 


-Im  x  (u)  = 


2  2 
(y  -  «<o  ) 

.  .2  2,  \ 

As  co(y  -  (o  )  ,  0  <  to  <  y 

,2  2, 

(y  -  a<o  ) 


(44a) 


(44b) 


(45a) 


=  ij  ttA  [ 6  ( <o  —  to  )  +  6(u>  +  u)_J  ]  ,  y  <  u)  <  “.(45b) 


Eqs.  (43)  and  (45)  satisfy  the  well-known  relation: 


2 

Im  xiu)  =  -  iTU)  aQ(w)  .  For  co  -*■  0,  Re  x(w)  =1  +  0  (to  )  ;  and  for 
to  ->  °°,  Re  x  (uj)  =  0(to-2)  aS  required.  In  the  high  frequency 
regime,  Im  x  (w)  agrees  with  Rajagopal.  For  a  =  1,  eqs.  43-45 
all  reduce  to  the  ideal  results  of  sec.  V. 

The  dynamic  structure  factor  can  be  obtained  by  the  relation 
it  S(to)  =  -  Im  x  (to)  .  Hence,  directly  from  (45)  we  have 
S  (to)  =  S(to)/x 


_  ,2  2, 

Ag  to  ( p  -  to  ) 

TT  (y2  -ctto2) 


0  <  to  <  y 


(46a) 


■=  ,jj  A  to  [5  (to  -  to  )  +  6  (to  +  to  )  ]  ,  y<to<».  (46b) 

P  P  P 

The  general  shape  of  the  dynamic  structure  factor  resembles  the 

three  dimensional  version  at  the  same  small  k,  consisting  of  a 

broad  single-particle  spectrum  and  a  sharp  collective  peak. 

However,  the  amplitude  of  the  two  dimensional  single-particle 

spectrum  is  "soft",  i.e.,  it  grows  with  a.  Hence,  there  is  an 

interesting  crossover  behavior  going  from  interacting  to  ideal, 
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wmch  is  unique  to  two  dimensions  at  small  k.  The  total  area 
under  the  single-particle  spectrum  remains  finite  at  all  values 
of  a  including  a  =  1.  The  dynamic  structure  factor  (46)  satisfies 
the  moment  sum  rules  to  all  orders  exactly  as  in  the  ideal  case. 


Intrinsic  Conductivit 


At  a  low  temperature  the  frequency-dependent  conductivity 
0^(01)  is  generally  understood  to  be  largely  due  to  impurities 
present  in  the  system.  The  electron-electron  interaction  con¬ 
tributes  to  it  hardly  at  all.  This  conventional  notion  ignores 
the  possible  existence  of  random  currents  arising  entirely  from 
density  fluctuations.  A  random  current  can  contribute  to  the 
frequency-dependent  conductivity  as  may  be  seen  from  the  Kubo 
conductivity  formula,  which  may  be  termed  an  intrinsic  con¬ 
ductivity.  It  appears  that  Kubo  had  speculated  upon  the 
existence  of  an  intrinsic  conductivity  some  years  ago.^ 

The  density  fluctuations  and  the  total  current  are 
connected  by  the  continuity  equation. 


p.  (t)  =  -  ik.j,  (t) 


(47) 


where  the  total  current  j^  includes  a  random  current  j^  .  Kubo 
showed  that  the  random  current  is  associated  with  the  random 
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force  F^ft)  such  that 


ik.j,  r(t)  =  F,  (t) 


(48) 


Since  our  system  contains  no  impurities,  the  Kubo  formula  assumes 


the  following  form: 

2 


ck(aj)  =  e 


..  “liOt  , 

dt  e  tj 


(jkr(t)  ,  jkr)/(jkr'jkr) 


(49) 


Hence,  using  (43)  one  can  calculate  the  intrinsic  conductivity 
present  in  a  pure,  homogeneous  system  if  the  random  force  is 
known.  This  contribution  must  presumably  be  added  to  the  con¬ 
ventional  mechanisms  to  obtain  the  total  frequency-dependent 
conductivity . 


Using  (13),  (15),  and  (31),  we  can  readily  evaluate  the 


random  current:  Suppressing  k-dependence, 

(jr(t),  jr)/(jr,jr)  =  b1(t)  =  2  J^ytj/yt  .  (5C 

2 

Hence,  writing  c  (u)  =  a  (u)/o,  where  a  =  2ire  p/my,  we  obtain. 


Re  a  U)  =  (y2  -  w2)  ^/y , 


0  <  w  <  y 


=  0,  y  <  w  <  °°. 


-Im  a  (s)  =u/y, 


0  <  a:  <  y 


=y/(u)  +  (u)2  -  y2)  1)  ,  y  <  u  <  oo.  (i 

The  real  part  of  the  conductivity  satisfies  the  conductivity 
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sum  rule: 


d\'  Se^liii)  =  p/2m. 


*  (53) 


We  note  that  the  intrinsic  conductivity  exists  in  low 
frequencies  only,  abruptly  vanishing  beyond  a  critical  frequency 
a  =  u  =  y  (see  Re  o,  (u) ) .  Detecting  the  intrinsic  conductivity 

C  K 
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raises  an  interesting  possibility.  In  an  impure  system,  e.g., 
metals,  in  addition  to  the  intrinsic  conductivity  there  is  the 
contribution  from  impurity  scattering.  Ordinarily  the  latter  is 
given  by  the  Drude  formula  which  is  a  high  frequency  approximation. 
Nonsmooth  behavior  in  the  vicinity  of  u>c  in  an  impure  system  would 
be  an  indication  of  merging  of  the  two  different  physical  mechanisms 


VIII.  Dynamic  RPA 


It  is  well  known  that  the  response  function  may  be  put  in 
the  form 

XkU,)  =  Xk(°W(l  +  Ak<u>)  xk  <0>  (o>) )  (54) 

where  Xk^°^  (<*>)  is  the  response  function  due  to  Hq  and  A^(aj)  is 

some  function  of  V  (see  (1)).  Various  dynamic  RPA  theories  are 
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equivalent  to  taking  A^(w)  =  A ( 0)  =  -  For  example,  A^  = 

gives  the  simple  RPA  and  =  v^d  -  G^)  gives  the  generalized 
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RPA,  where  is  a  local  field  term.  For  k«l,  *  O.  Hence, 
at  long  wavelengths  the  two  dynamic  RPA  theories  are  essentially 
the  same . 

Taking  A.  (gj)  to  be  frequency- independent  is  clearly  an 

A 

approximation,  the  validity  of  which  may  very  well  depend  on 
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the  static  regime.  Several  people  have  obtained  asymptotic 

conditions  in  three  dimensions  for  A^tu)  or  equivalently  for 
G^fw),  where  Ak(u)  =  v^  (1  -  G^fw)),  sometimes  referred  to  as 
a  dynamic  local  field  term.  There  are  otherwise  no  exact 
general  expressions  known.  Using  the  method  of  recurrence 
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relations  we  have  recently  derived  an  expression  for  Ajt(z=itu): 
Suppressing  k  dependence 

A  ( z)  =  [x*1-  (x(0))'1]+  1  (bi<z>>“1  "  (b1(o)  (z) )  ~1] 

=  A  +  )  (z)  .  (55) 

The  first  bracketed  term.  A  is  z-independent .  In  general, 
b^(z)  t*  b^°^(z)  owing  to  the  fact  that  ,  v  >  2.  A 

change  in  the  memory  function  when  interaction  is  turned  on  makes 
the  local  field  term  frequency-dependent. 


In  the  two  dimensional  electron  gas  at  long  wavelengths, 

we  have  shown  that  b^(z)  =  (z)  ,  since  ^v=  »  v  >  2 

(see  sec.  VI) .  That  is,  when  density  fluctuations  are  confined 

to  the  vicinity  of  the  Fermi  surface,  the  memory  function  is 

unaffected  by  the  electron-electron  interaction.  Hence,  in  this 

static  regime  i.e.,  k<<l,  A^(ai)  =  and  the  dynamic  RPA  theorie 
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are  valid.  This  accounts  for  the  RPA  work  of  Holas  et  al. 


IX. 


Discussion 


The  method  of  recurrence  relations  treats  the  time  evolution 
of  dynamical  variables  as  a  geometric  problem  in  Hilbert  space. 

In  this  approach,  the  recurrants  are  basic  elements  of  a  realized 
Hilbert  space  whose  structure  contains  the  possible  time  evolu¬ 
tion.  The  recurrence  relations  RRI  and  II  are  merely  mathe¬ 
matical  expressions  of  this  geometric  picture. 

For  the  two  dimensional  electron  gas  at  long  wavelengths, 
the  recurrants  assume  a  very  simple  form  making  the  structure 
of  the  Hilbert  space  very  elementary.  As  a  result,  we  were 
able  to  deduce  an  exact  solution  for  the  time  evolution  of 
density  fluctuations  occurring  near  the  Fermi  surface.  Using 
this  solution,  we  have  also  obtained  other  physical  quantities 
via  linear  response  relations.  In  some  other  static  regimes 
e.g.,  k  21  kp,  the  recurrants  should  have  some  other  forms. 

Hence,  we  expect  the  time  evolution  of  density  fluctuations  to 
be  different. 

Even  when  confined  to  the  vicinity  of  the  Fermi  surface, 
the  density  fluctuations  show  persistent  memory  effects.  There 
thus  exists  a  random  current  and  there  is  an  intrinsic  con¬ 
ductivity.  Remarkably,  the  memory  is  shaped  by  quantum 
fluctuations  and  not  by  the  electron-electron  interaction. 

This  has  an  unexpected  consequence  that  the  dynamic  RPA  theories 
are  exactly  valid  in  this  static  regime. 

The  time  evolution  of  density  fluctuations  may  be  regarded 
as  a  Fermi  sphere  undergoing  a  restoring  motion  due  to  a  slight 
initial  displacement  in  momentum  space.  In  the  two  dimensional 


electron  gas  there  are  an  infinite  degrees  of  freedom  for  such 
a  relaxation,  giving  rise  to  an  infinite  dimensions  for  the 
corresponding  Hilbert  space.  The  resulting  density  fluctuations 
are  electron-hole  scattering  processes  on  which  a  collective 
process  is  superimposed.  These  processes  are  dimension-sensitive 
In  one  dimensional  electron  gas,  for  example,  the  restoring 
motion  is  one  dimensional  and  the  single-particle  and  collective 
processes  are  bound  together.  For  this  system,  the  Hilbert 
space  is  also  finite-dimensional  making  the  relaxation  always 


oscillato 
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Appendix  A.  Relaxation  function 


We  shall  now  evaluate  the  integration  given  below: 


ao(t)  •  KT  /  dz  e 


zt  (Z2+U2)  "  -  (1-s) z 


2  2, 

z  +y  /a 


(Al) 


_  u 

Since  a<l,  there  are  isolated  poles  at  z  =  ±iya  =  ±i(jp,  which 


are  outside  the  cut  bounded  by  z  =  ±iu  on  the  imaginary  axis.  Hence 
the  contributions  from  the  isolated  poles  are  straightforward  to 


evaluate.  The  residue  at  z  =  ±icj  is: 

P 


s/2  a'  [/I  -  (y/w  >*  +  (1-s)]  e 

Lr 


±iaj_t 


=  1/1=5  -  (l-a) ] e 


±iw  t 

P  .  (A2) 


The  integration  about  contour  c.^  containing  the  cut  is  as  follows: 


s/2~iy’/  dz  e 


7  7  ^ 

zt  (y  +  z^)  * 


1  +  a  z2/y 2 


=  s/yt*/  dv  /y2  -  y2  [1  +  ay2/y 2  -f  a2y4/y 4  +  .  . .  1  e^"yt 


/_  r  ,  ?  ri  a  ,9  \2,gl  /9  v  4 

sA;,-/  dy  /y  -y  [1  -  ^ ^  ' 

—y  y  y 


.  ]e 


iyt 


^  «  2  ^  4 

j,  r  i  Q  /  °  \  2  oc  #  9  \ 

[1  “  2(St'  4  St 

y  y 


.  •  ]  Jj,  (yt)  /yt 


(A3) 


where  we  have  used  an  integral  representation  of  the  Bessel 


function 


Combining  the  above  results  we  have  with  the  definitions  (42a ,b) 

(A4) 


a  (t)  =  A  £(-a)V  (-1— r-)  2v  J.(yt)/yt  +  A  cos  w  t. 
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Appendix  B.  Proof  that  aQ(t=0) 

We  shall  consider  t-*0  in  our  solution,  for  a  (t) 

o 

Appendix  A: 


a  (t)  *  A  l  {-a)v  (|— r)2j  J-,(yt)/ut  +  A  cos  ai, 
°  sv=0  1  p 

Using  the  expansion  formula,  we  have 

oo 

(z)/2  =  1  C2k  2  * 

1  k=0  ZK 

where 

C2k  =  ^)k/(22k+1  k!  (k+1)  1)  . 

Hence, 

‘I?’2"  JX<2)/2  !  =  2vI  e2u  5  «2v- 

Substituting  these  results  in  (Bl) ,  we  obtain 

ao<°>  *  As  52V  *  V 

n=0 

Now  one  can  show  that 


OO  \) 

v  Ct  2v !  _  171  f\  •  _■ , 

£  2-. > _ i  —  r  ( 1  ~  'i  “  a )  . 

v=0  2^  x  v!  (v+1)  !  a 


Using  A  =  1  -  vl  -  a  ana  A  =  2(/l  -  a  -  (1  -  a)) /a, 


found  in 


t.  (: 


(] 


(] 


(3 


(I 


we  find 


Appendix  C.  Proof  that  a  (t)  =  (t)  when  a  =  1 


We  shall  consider  a-1  in  our  solution  for  aQ(t)  given  below 

CD 

a  (t)  =  A  l  (-)V  z-)  V  J,  (y  t)  /y  t  +  A  cos  w  t.  (Cl) 

°  Sv-0  O  M  1  P  P 

For  a-1 ,  we  can  write  setting  Ag  =  1  and  A^  =  0 , 

a  It)  =  l  d.  (a-1)  (ut) 2k.  (C2) 

°  k=0 

Although  d^{a)  is  complicated,  it  becomes  simple  when  a  =  1: 

dk(l)  =  (-)k/(22k  <k!)2)  .  (C3) 

Hence, 

00  ,  ,k  ,  ,  ,2k 

lim  a  (t)  =  l  >Z>  -  =  j  (ut).  .  (C4) 

a-1  0  k=0  (kir  0 


We  recall  that  av  !  (t)  =  J  (ut) ,  comoleting  our  proof. 

o  o  * 


Appendix  D.  Frequency-Dependent  Relaxation  Function 


We  shall  first  calculate  the  frequency-dependent  relaxation 
function  due  to  single-particle  excitations  a® (u)  via  the  formula 

®  '  •  || 

aS(io)  =  1/277-  /  dt  elujt*eitl  aS(t)  (Dl) 

°  _»  ° 

where  e-*-0.  Using  (Bl)  ,  we  have 

_  S  /  ,  \  _  »  T  /  2»v  (2v)  fn'iS 


a“('a!)  =  A  /2tt •  £(-a/yz)v  <D 
°  n=0 


where 


(2 v)  _  .3  ,2v  r  -i-t  -  e  t| 


-  (fe)"  /  dt  e- 


(tit)  /lit. 


,  (o)  _  .  2  2.  k  _ _  ^ 

=  0  /  y<u<®  .  (D4 ) 

Now  integrating  by  parts,  we  get  <?^.2v^  *  (iu;)  2v  <{> .  Hence,  using 
these  results  we  obtain  directly 

aS  (yj)  =  —  A  (u2  -  u>2)  ?V(1  -  ota2/y2)  ,  0<o»<y 

O  77  S 

=  0  U<(i}<“.  (D5) 

The  relaxation  function  due  to  the  plasmon  excitations  a^(u)  is 
straightforward: 

aP(x)  =  1/2  A  (5(u;-u  )  +  5(w  +  w  )).  (D6) 

o  p  p  p 

Since  u  =  u/a2,  where  a<l,  the  plasmon  contributes  only  for 


Appendix  E.  The  Response  Function  x(^) 


We  shall  use  the  well-known  relation  given  below  to  calculate 
the  response  function  from  the  relaxation  function, 

CO 

X(a)  »  1  -  iu  /  dt  e"1U)t  “  £t  aQ(t)  (El) 

o  '' 

where  e-*0 .  From  our  solution,  we  can  write: 

ao(t)  5  VI(t)  +  A0N(t)*  (E2) 

One  can  immediately  show  that 

U(io)  =  ia)/(y2  -  u2)  +  ri(5(w  -w)  +  5  (u  +w)  )  ,  (E3) 

P  2  r  p 

•  u 

where  w  =  u/a  .  To  calculate  M(w)  ,  we  first  define 
P 

r  U)  =  /  dt  e'lUJt  “  £t  (It-t)2"’  J,  (ut)/ut  (E4) 

^  g  .  OUt  1 

Then , 

CO 

M(w)  =  I  (-a)V  F-U).  (E5) 

n-0  z 

Now  integrating  by  parts,  we  find  a  recurrence  relation  for  T's: 


=  -iw  a  -  u 

;  2v-2 


2v-2 


,  v>l  , 


where  =  (|— ■ r)2^  J.(yt)/ut.  The  above  recurrence  relation  can 

^  t1  0  t  -L 

be  used  to  sum  (E5) .  By  substituting  (E6)  in  (E5) ,  we  get 


2  2  \) 

M  =  7  +  a  a;  M/u  +  iaw  J  {-a)  g~ 

°  v=0 


M  ■  : - T-T  ‘ro  +  I  g2v'- 

1  -  Otiii  /  y  v  =  0 


The  sum  has  been  already  carried  out  (see  Appendix  B) .  Also  T 
can  be  evaluated: 


F  =  [  ( u  2  -  w2)  ^  -  iai]/u^  ,  0<u><y 

o 


p<U)<«  . 


(] 


Together 

M=  (^-au2)"1  (/u2  -  ^  -  i'jj/1  -  a)  ,  0<uj<u 

=  i  (u2  -  aw2) _1  (/jj-u  -oj/1  -  a)  ,  U<w<*>.  (] 

Combining  (E2)  and  (E10) ,  we  get 


Re  V  =  i  +  w  A„  Im  M  +  u  A  Im  N 
A  s  P 

-  1  +  A  c,2/(l  -~u)  /  (u2  -  au2)  ‘  0<U)<P 

s 

=  1  +  As  u)2IvT^  +  A  -  u2/A)/(y2  -  aw2)  ,  y<w<~  C 
where  we  have  used  A^  —  2A^  >'1  -  a/a.  Similarly/ 


jj  A  Re  H  T  i  A_  Re  N 
s  P 

,  2  2  ...  2  2. 
t\  CJ  v  U  ~  u)  /  (  -i  “Ct-0  )  , 

S 

0<w<y 

Utt  Ag  iu  ( 5  (i->  —  ^p)  +  5  (u  +  Up)  )  / 

y  <U)<  »  . 

( 

Observe  that  our  results  satisfy  the  relationship  Im  x  -  -uua 
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